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We review the current status of small x resummation of evolution of parton distributions and of deep-inelastic 
coefficient functions. We show that the resummed perturbative expansion is stable, robust upon different treat- 
ments of subleading terms, and that it matches smoothly to the unresummed perturbative expansions, with 
corrections which are of the same order as the typical NNLO ones in the HERA kinematic region. We discuss 
different approaches to small x resummation: we show that the ambiguities in the resummation procedure are 
small, provided all parametrically enhanced terms are included in the resummation and properly matched. 
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1. The need for small x resummation 

The so-called small x regime of QCD is the 
kinematical region in which hard scattering pro- 
cesses happen at a center-of-mass energy which 
is much larger than the characteristic hard scale 
of the process. An understanding of strong in- 
teractions in this region is therefore necessary to 
do physics at high-energy colliders. In this sense, 
HERA was the first small x machine, and LHC is 
going to be even more of a small x accelerator. 

In deep-inelastic lepton-hadron interactions, 
the scale is set by the virtuality of the photon 
Q 2 , and x — ^— = + 0(x)), where s is 

the center-of-mass energy of the virtual photon- 
hadron collision, and the distribution of partons 
which carry a fraction x of the incoming nucleon 
energy is probed. In hadron-hadron interactions, 
the scale is set by the invariant mass M 2 of the 
final state, and x — X1X2, with Xi = M=e ±v , s the 
center-of-mass energy and y the rapidity of the 
hadron-hadron collision. Here, the distribution 
of partons which carry fractions X{ of the two in- 
coming nucleon energies are probed. This means 



that typical x values probed at the LHC in the 
central rapidity region are almost two orders of 
magnitude smaller than x values probed at HERA 
at the same scale. Hence, small x corrections 
start being relevant even for a final state with 
a characteristic electroweak scale M ~ 100 GcV 
(see Fig. [[]). 

As is by now well known, perturbative correc- 
tions become large at small x. Due to the ac- 
cidental vanishing of some coefficients, the lead- 
ing large corrections cannot be seen in NLO and 
NNLO splitting functions; however, the first sub- 
leading correction can already be seen in the 
NNLO splitting functions which have been com- 
puted recently, as well as in NNLO coefficient 
functions: they are large enough to make recent 
NNLO parton fits unstable at small x (see Fig. [5]). 

This suggests dramatic effects from yet higher 
orders, so the success of NLO perturbation theory 
at HERA, as demonstrated by the scaling laws 
it predicts [314] . has been for a long time very 
hard to explain. In the last several years this sit- 
uation has been clarified, mostly thanks to the 
effort of two groups (ABF [S]-[1Q] and CCSS p]- 
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Figure 1. Comparison of the HERA and LHC 
kinematical regions (from Ref. pQ). 



[16]), which have presented a full resummation of 
evolution equations in the gluon sector, thereby 
showing that, once all relevant large terms are 
included, the effect of the resummation of terms 
which are enhanced at small x is perceptible but 
moderate — comparable in size to typical NNLO 
fixed order GLAP corrections in the HERA re- 
gion. A detailed comparison of the ABF and 
CCSS approach in the pure gluon (n/ = 0) case 
was presented in Ref. [1 , where excellent agree- 
ment was found. 

This approach has now been generalized by 
ABF [17j to a full resummation including quarks, 
and including the resummation of deep-inelastic 
coefficient functions, so that resummed expres- 
sions for deep-inelastic structure functions can 
be obtained. Progress towards the inclusion 
of quarks has also been made by the CCSS 
group [18], though full results are not yet avail- 
able. Meanwhile, an alternative somewhat simpli- 
fied approach has also been developed [15] . and 
used to perform a fit to deep-inelastic scatter- 
ing data [20 . In this approach, the factoriza- 
tion scheme is not defined in a fully consistent 
way at the resummed level, and also some con- 
tributions to the resummation arc cither not in- 
cluded, or treated by means of truncated expan- 
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Figure 2. Comparison of the LO, NLO and NNLO 
gluon distributions in the MSTW08 parton fit 
(from Ref. 0). 



sions. A comparison of resummed predictions for 
deep-inelastic structure functions obtained in the 
TW and ABF schemes for resummation is pre- 
sented in Ref. [2T] , 

Here, we shall review the main ingredients that 
go into small x resummation, with the aim of un- 
derstanding the impact of various contribution on 
the final result and its stability upon higher or- 
der perturbative corrections and upon the inclu- 
sion of various subleading contributions: we shall 
take the ABF approach as a baseline, and discuss 
the impact of alternative options. We shall show 
that small x resummation is very constrained by 
various requirements, which include momentum 
conservation, the inclusion of collincar contribu- 
tion and matching to GLAP evolution, and con- 
sistency with the renormalization group. Once 
all these requirements are met, further sublead- 
ing ambiguities are quite small; however, if some 
of the corresponding terms are missed out, effects 
are not negligible as we shall see. The resummed 
corrections thus obtained are perturbatively sta- 
ble, as demonstrated by the fact that renormal- 
ization and factorization scale dependence are 
moderate, and decrease with increasing pertur- 
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bative order as they ought to. The typical effect 
of the resummation in the HERA and LHC re- 
gions is comparable to that of NNLO corrections, 
but with the opposite sign. 

In the next section, we shall review the ingre- 
dients which are necessary in order to perform 
the resummation in the gluon sector. In the sub- 
sequent section, we shall summarize the generic 
features of the resummed results. In the last 
section, we shall discuss how quarks and deep- 
inelastic coefficient functions may be included in 
the resummation, and discuss resummed results 
for deep-inelastic physical observables. 

2. The three ingredients of stable resum- 
mation 

In this section, we discuss the resummation of 
evolution equations when rif = 0. In this case, 
there is a single parton distribution, the gluon 
distribution G(£,t), with £ = In A, t = ha^. It 
is convenient to define the Mellin transforms 

/>oo 

G(N,t) = / d£e- Ni G(£,t) 
Jo 

/oc 
dte- Mt G(Z,t) (1) 
-oc 

2.1. Double leading expansion 
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Figure 3. Double leading expansion of the GLAP 
anomalous dimension 7 (left) and the BFKL ker- 
nel x (right). 

The gluon distribution G(N, t) can be ex- 
pressed in terms of the gluon distribution at 



t = < by solving the GLAP equation 

j t G(N,t)= 1 (N,a s )G(N,t) (2) 

which at the N fc LO sums all terms of order 
a™t n ~ k , to all orders in a s . The first step of re- 
summation consists of including, to the N fc L log 
level, all contributions to the anomalous dimen- 
sion 7(iV, a s ) of order a™N~t n ~ k \ to all orders 
in a s , since they correspond to contributions of 
order a™ ln"- fc ± to G(£,t). 

This inclusion is straightforward at the fixed 
coupling level, thanks to the fact that the gluon 
distribution G(£, M) can be expressed in terms 
of the gluon distribution at £ = £0 by solving the 
BFKL equation 

±G(£,M) = X (M,a a )G(Z,M), (3) 

whose kernel xiMi a s) IS simply the inverse 
function of the GLAP anomalous dimension 
l(N,a 8 ) [US]: 

xh(N,a s ),a s ) = N. (4) 

The duality equation maps the perturba- 
tive expansion of 7 in powers of a s at fixed N 
in the expansion of x i n powers of a s at fixed M 
and conversely. One can thus construct a double 
leading expansion [22] (see Fig. [3]) which includes 
in x an terms up to a given order in the expan- 
sion in powers of a s both at fixed M and at fixed 
jj. Its dual 7 can be shown to include terms up 
to the same order in the expansion in powers of 
a both at fixed a s and at fixed jf. 

Using either the double-leading x or the double 
leading 7 in the BFKL or GLAP equation respec- 
tively leads to a solution which includes all terms 
which are logarithmically enhanced either in i 
or in Q 2 to the given order [B]. The result (see 
Fig. @]) is close to the GLAP one when M — ► 0, 
and close to the BFKL one when N — > 0. Because 
the perturbative expansion of the BFKL kernel is 
very poorly behaved, this resummed result has 
poor perturbative stability as N — > 0. 

2.2. Exchange symmetry 

The perturbative instability of the kernel as 
N — > can be cured by observing that the BFKL 
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Figure 4. The BFKL kernel and its dual GLAP 
anomalous dimension computed at LO and NLO 
in the BFKL expansion, the GLAP expansion and 
the double-leading expansion. 

kernel must be symmetric upon the interchange 
M — > 1 — M, due to the symmetry of the three- 
gluon vertex upon the interchange of the radiated 
and radiating gluon [13j . This symmetry, which 
is manifest in the LO BFKL kernel, is however 
broken beyond the LO of the BFKL expansion 
by running coupling correction and by the choice 
of DIS kinematics [11I13| . Nevertheless, the sym- 
metry breaking terms can be computed exactly. 
Hence, one can symmetrize the double-leading 
expansion by undoing the symmetry breaking 
terms (by changing kinematics and argument of 
the running coupling) , then symmetrizing the re- 
sults, and finally restoring the original symmetry- 
breaking kinematics and choice of argument of 

The result turns out to be surprisingly stable 
because of two features: (a) the anomalous di- 
mension must satisfy the momentum conservation 
constraint j(l,a s ) = 0, which by duality implies 
x(0, a s ) = 1; the anomalous dimension j(N, a s ) 
decreases monotonically as N increases as a con- 
sequence of the fact that a gluon looses momen- 
tum when radiating. Combining these two, it 
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Figure 5. The LO and NLO resummed sym- 
metrized double-leading kernels compared to the 
LO and NLO kernels in the BFKL expansion and 
the NLO GLAP kernel. CCSS denotes the corre- 
sponding result of Ref. 15J (from Ref. [1]) 

follows that x always has a minimum, because 
it must go through the value x = 1 at a pair 
of values of M which are symmetric about the 
minimum, located at M = \. The minimum 
is preserved even when the symmetry is broken, 
in which case the two "momentum conservation" 
points at which x — 1 get shifted to M = and 
M — 2. One can further show that when the 
symmetry breaking is removed, the kernel is an 
entire function in the M plane; after restoring the 
symmetry breaking it remains free of singularities 
for ReM > -1 [10]. 

The perturbative instability of the BFKL ker- 
nel is thus completely removed: see Fig. [5j In 
this Figure we also compare results with those 
obtained by CCSS through a procedure which 
is rather different, but shares the following fea- 
tures: (a) all logarithmically enhanced terms in 
Q 2 and i are included up to NLO (b) the under- 
lying symmetry is implemented up to sublead- 
ing terms. The extreme similarity of the results 
demonstrates the stability of the procedure. 

2.3. Running coupling 

The double-leading expansion upon which the 
resummation has been based so far sums up all 
terms which are large when a s is small, but 
a s {; — a s In — ~ 1: namely, a contribution of the 



Can we trust small x resummation? 



5 



form F(a s £) to the splitting function is consid- 
ered to be of order a°. Thus, if the computation 
is performed at N fc L order, relative corrections 
are of order a^ +1 when a s — > while either x or 
a s £ are kept fixed. However, this does not guar- 
antee that corrections remain small if £ — > oo (i.e. 
x — > 0) with a s fixed. It is easy to see [5] that 
any correction which changes the leading (right- 
most in the N plane) small N singularity of the 
anomalous dimension leads to a contribution to 
the splitting function which blows up as £ — * oo 
in comparison to the lower order, because such a 
contribution changes the asymptotic x — » be- 
haviour of the splitting function, and thus of the 
parton distribution. 

The leading singularity of the anomalous di- 
mension is a simple pole at TV = in the LO (and 
NLO) GLAP case. After double leading resum- 
mation it becomes a square- root branch cut: the 
inverse dual Eq. (HJ of the quadratic behaviour 
of the kernel near its minimum. As discussed 
in Sect. 12. 2i the presence of a minimum of the 
kernel is a generic feature which follows from its 
symmetry properties and momentum conserva- 
tion. The intercept of the kernel however changes 
order by order, and it is only its all-order posi- 
tion which determines the asymptotic small x be- 
haviour. Hence, higher order corrections to the 
position of the minimum are asymptotically large 
as £ — > oo: this suggests that the all-order lo- 
cation of the minimum of the kernel should be 
treated as a non-perturbative parameter, to be 
fitted to the data |6|7j . 

However, running coupling corrections change 
this state of affairs. Because the running of the 
coupling is subleading in In i , running coupling 
corrections only affect the duality equation ([4} 
through a finite number of terms: at the N fc L 
logarithmic level, the duality relation is corrected 
by a term of the form 

1 { £ = M)"fff (f ) , (5) 

where the function fj^ (jf) can be calculated at 
any given order using suitable operator methods 
in terms of the N fc LO BFKL kernel [23] . 

Now, it turns out that whenever the fixed cou- 
pling kernel has a minimum, the running cou- 



pling correction 7k ' is also asymptotically large 

as £ — > 00: in fact f$ Eq. §5§ grows as in 
comparison to the splitting function computed at 
LO in the double-leading expansion. Hence, in 
order to determine the correct £ — * 00 limit we 
must resum all these terms. This resummation 
can be performed exactly for the series of terms 
which grow fastest as £ — > 00 [8]: the result can 
be expressed in terms of Airy functions for a ker- 
nel which is linear in a s , and Bateman functions 
for a kernel with a generic dependence of a s [10] . 
The (asymptotic, divergent) expansion of these 
functions in powers of a s at fixed % gives back, 
to each order in a s , the contribution which upon 
inverse Mellin transform grows fastest as £ — > 00 
to the terms 7^ Eq. ([5]). 
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Figure 6. Quadratic approximation to the NLO 
resummed kernel of Fig. [5] and its Bateman 
running-coupling resummation. 

It can then be seen that the resummation 
of running coupling corrections changes the na- 
ture of the leading singularity: the fixed-coupling 
square-root branch cut is turned back into a sim- 
ple pole. This is shown in Fig. [6] where we dis- 
play the quadratic approximation to the double- 
leading NLO kernel of Fig. [5] and its Bateman re- 
summation, i.e. the anomalous dimension which 
is obtained from it when running coupling correc- 
tions to Eq. ([5]) it are included to all orders (com- 
puted using the dependence on a s of the NLO 
resummed result of Fig. [5]). 
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Hence, after running coupling resummation, 
the minimum of the kernel no longer provides the 
leading small x singularity, which is instead given 
by the pole in the Bateman anomalous dimension. 
The location and residue of the Bateman pole are 
fully determined by the intercept and curvature 
of the minimum of the original kernel, and their 
dependence on a s . 

3. General features of resummed results 
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Figure 7. The resummed splitting function with 
?i/ = and a s = 0.12 with fixed coupling (top) 
and running coupling (bottom). 

The crucial feature of the resummation proce- 
dure summarized in the previous section is that 
at each step the contributions which are included 



in the anomalous dimension on top of its stan- 
dard GLAP fixed-order expansion in powers of 
a s at fixed N can be expanded out perturba- 
tively, so that it is possible to obtain a fully re- 
summed expression by simply combining all con- 
tributions, and subtracting the double counting 
terms. The procedure can be performed order by 
order in perturbation theory, by starting with the 
double- leading expansion of Fig. [3] to any given 
order, and then improving it as discussed in the 
previous section. 

The resummed anomalous dimension has then 
schematically the form 

lZ C NLoMt),N) = T^ffoW^AO 
+ 1 B (a s (t),N)~ 1 ?(a s (t),N) 

-Tf.M*),^-<oM*),#) (6) 

+7match (««(*), N) + Jmam(ot a (t),N), 

where 

• ^NLo( as (t)' is the fixed-coupling re- 
summed anomalous dimension displayed in 
Fig. El obtained by duality Eq. (j4|) from 
the kernel which in turn is found by sym- 
metrization (Sect. |2~2")) of the NLO double- 
leading kernel (Figs. EH}; 

• 7 B (a s (t) , N) is the Bateman anomalous di- 
mension Fig. [SI 

. 7 ? M*), AO, lf s (a s (t),N) j? St0 (a s (t),N) 
are double counting subtractions between 
the previous two, namely the contributions 
to the LO and NLO terms Eq. © 
which grow fastest as £ — ► oo; 

• 7mom subtracts subleading terms which 
would otherwise spoil momentum conserva- 
tion; 

• 7match subtracts any contribution which de- 
viates from NLO GLAP and at large N 
(which corresponds to large x) doesn't drop 
at least as jj. 

Note that the last two contributions are formally 
subleading: they are included in order to im- 
prove the matching to the GLAP anomalous di- 
mension, in that they remove from the resumma- 
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tion subleading contributions which may be non- 
negligible in the large x region where the resum- 
mation is not supposed to have any effect. 

In Figure [7] we display the splitting functions 
obtained from Mellin inversion of the resummed 
anomalous dimension Eq. |J7J at the fixed cou- 
pling level and at the running coupling level — 
i.e. respectively without and with the inclusion 
of the Bateman contribution j B (a s (t), N) and its 
associate double counting subtractions. The re- 
sult is also compared to the CCSS result Ref. [15] 
(from Ref. pQ). The resummed expansion is seen 
to be stable (the LO and NLO results are close), 
all the more so at the running coupling level. The 
resummed result matches smoothly to the GLAP 
result in the large x < 0.1 region, but at small x 
it is free of the instability which the GLAP ex- 
pansion shows already at NNLO. 

The comparison between results obtained with 
the ABF method discussed here and those by 
CCSS [15] is illuminating in various respects. The 
CCSS approach also includes the three ingredi- 
ents discussed in the previous section — double- 
leading resummation, symmetrization of the ker- 
nel, and running coupling corrections — but the 
implementation is rather different. In particular, 
the resummation of running coupling corrections 
is obtained by numerical solution of the running- 
coupling BFKL equation J3]) (see Ref. []]). The 
closeness of results obtained by CCSS and ABF 
at the fixed coupling level reflects the closeness 
of the kernels Fig. [5] The fact that CCSS and 
ABF results are even closer at the running cou- 
pling levels follows from the softening of resum- 
mation effects due to asymptotic freedom. Also, 
the fact that exact numerical resolution of the 
running-coupling BFKL equation (in the CCSS 
approach) followed by a numerical extraction of 
the associate anomalous dimension, and the re- 
summation of the leading running coupling cor- 
rections Eq. (0 in terms of a Bateman function 
(in the ABF approach) lead to a result which 
manifestly coincides for all x 0.03 supports the 
accuracy of both procedures. 

It is interesting to observe however that the 
CCSS resummed result shows a significant devi- 
ation from the unresummed GLAP result even 
for x iS 0.1, which is not seen in the ABF re- 



sult. Because this effect is only present at the 
running coupling level, it is likely to be due to 
the fact that running coupling terms lead to con- 
tributions which survive in the large x region, 
where these terms are formally subleading (recall 
that the running coupling contributions were se- 
lected on the basis of their behaviour as £ — > oo 
i.e. x —y 0). In the ABF approach (but not 
in the CCSS approach) these contributions are 
subtracted through the term 7 ma tch (c*s {t) , N) in 
Eq. ([7]). We have found this subtraction to be 
necessary in order for the resummed results not 
to differ from the fixed-order ones in the x ^ 0.1 
region. 



. lead sing / 

/ 

0.5 — BFKL LO / 

/ sym NLO 




-0.3 



W 1 2 3 4 5 6 

Figure 8. The leading small x singularity vs. a s 
(top) and the slope of the splitting function vs. 
£ = In i (bottom, from Ref. [25]). 

More detailed insight into the generic features 
of the resummation can be obtained by study- 
ing its small x behaviour. As already mentioned, 
this is largely determined by the position of the 
leading (rightmost) singularity of the anomalous 
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dimension in the N plane. For the LO and NLO 
GLAP anomalous dimension this is located at 
N = 0, while at the resummed level it is plot- 
ted as a function of a s in Fig. [5J the large and 
perturbatively unstable corrections at the BFKL 
level (compare Fig. [4]) are turned into more mod- 
erate and stable correction (compare Fig. [5]) by 
the symmetrization, and further reduced and sta- 
bilized by the running of the coupling. 

The relative importance of various contribu- 
tions to the resummation is elucidated by com- 
paring in Fig. [8] the slope of the resummed 
running-coupling splitting function (Fig. [7|) with 
respect to £ to that of the NLO GLAP result, 
and that of the splitting function obtained by in- 
verse Mellin transform of the Bateman anomalous 
dimension j B (a s (t), N) in Eq. |[7J). It is clear 
that for x <> 0.1 the slope (and thus, up to a 
constant, the splitting function) is determined by 
the GLAP result, and for x Sa 10~ 2 by the Bate- 
man result [25] . Note that this is despite the fact 
that it is only for tiny values of x Si 1CP 15 that 
the Bateman asymptotic slope reduces (to within 
about 10%) to its asymptotic value No (the loca- 
tion of the leading singularity) [2"5] . 

Hence the qualitative features of the resumma- 
tion are essentially determined by matching to the 
GLAP result the Bateman anomalous dimension, 
which in turn is fully determined by resummation 
of the running coupling terms Eq. (0, and thus 
by the intercept and shape of the minimum of 
the symmetrized kernel displayed in Fig. [5] [lOj . 
The insensitivity of these to the details of the re- 
summation procedure explains the stability of the 
results. The dominant feature of the result in the 
HERA region 10~ 2 ^ x Z 10 -4 is a dip in the 
splitting function which results from this match- 
ing (see Ref. [16] for an alternative discussion of 
this feature of resummed results). 

A resummed gluon splitting function was also 
presented in Refs. [19120] . The agreement is rea- 
sonable at the qualitative level, but the resummed 
splitting function appears to display a stronger 
rise at small x and a somewhat more pronounced 
dip at intermediate x. This may be due to the fact 
that in Ref. |19|20j no symmetrization of the ker- 
nel is performed: this, as discussed above, leads 
to resummed results which tend to be unstable 



upon pcrturbative corrections as x — » 0. 

4. Quarks and phenomenology 

Going from a resummation of the evolution 
equations at Nf — to fully resummed physical 
observables requires two ingredients: the inclu- 
sion of the quark sector in the resummation, and 
resummed partonic cross sections. 

4.1. The resummed splitting function ma- 
trix 

An extension of small x BFKL-like evolution 
equations to the quark sector has been attempted 
in Ref. [18], but results in a fully consistent fac- 
torization scheme are still not available. How- 
ever, it turns out that this construction is not 
necessary for the determination of the matrix of 
resummed splitting functions. This is due to the 
fact that the resummation only affects one of the 
two eigenvectors of the singlet anomalous dimen- 
sion matrix. Therefore, in order to obtain coupled 
resummed evolution equations for singlet quarks 
and gluons it is sufficient to fix the factorization 
scheme at the resummed level [17] , 

Indeed, calling 7+ and 7" the two eigenvectors 
of the anomalous dimension matrix, if only con- 
tributions to 7 which are singular as N — > are 
included, 7~ = and thus 

up to nonsingular terms. Furthermore, in the MS 
and related schemes, 7" = jf^l" (up to nonsin- 
gular terms) [27] . Hence, combining the determi- 
nation of 7 + with the knowledge of 7 9S in the MS 
scheme, which was computed at the resummed 
level in Ref. [27], the anomalous dimension ma- 
trix is fully determined. 

The actual computation of the full splitting 
function matrix with nj ^ entails some techni- 
cal complications, which have only recently been 
solved in Ref. [17]. Firstly, when n/ 7^ 0, the 
eigenvectors of the anomalous dimension matrix 
develop an unphysical singularity at the value of 
N where the two eigenvectors cross. This singu- 
larity cancels in the solution to evolution equa- 
tions, and the cancellation must be enforced ex- 
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Figure 9. The splitting function matrix with 
rif = 4, a s = 0.2. The curves correspond to 
LO (dotted), NLO (dashed), NNLO (dot-dash), 
resummed (solid). The two different resummed 
curves (in the gluon sector) correspond to the 
MS (steeper at small x) and QoMS factorization 
schemes. 



actly throughout the double-leading resumma- 
tion and symmctrization: if it were spoiled by 
subleading terms, this would lead to a spurious 
small x rise of parton distributions. Furthermore, 
the quark-sector anomalous dimension 7 gff was 
determined in Ref. [57] in the MS scheme, while 
the double-leading resummation is most natu- 
rally performed in the Qo scheme [29128] : in the 
MS scheme the running coupling terms Eq. ([5]) are 
factored order by order in the coefficient function. 
Because their resummation determines the lead- 
ing small x behaviour, as discussed in Sect. 12. 3[ it 
is more convenient to perform the resummation 
in a scheme where they are included in the split- 
ting function. Whereas the scheme change from 
the Qo scheme to MS in the pure-gluon case was 
worked out previously [29)28] . its construction in 
the presence of quarks is nontrivial [22|17j . 

The resummed matrix of splitting functions 
in the MS and QoMS factorization schemes is 
compared to the unresummed result in Fig. [9] 
Whereas a detailed comparison with the CCSS 
splitting function matrix [18] is not possible, 
because CCSS results are given in a scheme 
which differs from the standard MS scheme by 



an amount which is undetermined beyond fixed 
NLO, a qualitative comparison shows reasonable 
agreement. 

Quark sector splitting functions have also been 
given in Ref. [20] (see also Ref. [21]). Their agree- 
ment with those of Fig. [9] is not so good: the P qg 
splitting function shows a much stronger small x 
rise, and a sizable deviation from the NLO re- 
sult at large x ^ 0.1. The latter feature can be 
understood as a consequence of the fact that in 
Refs. [19120] no matching to large x is performed: 
contributions from the small x resummation in 
the large x region are not subtracted. The for- 
mer feature is likely to be due to the fact that 
results in Ref. [20] are not determined in a fully 
consistent factorization scheme. In fact, the re- 
summation is performed in the Qo scheme, but it 
is then combined with MS (or DIS) scheme coef- 
ficient functions: in the TW approach, the issue 
of the scheme transformation from Qo to MS is 
still unresolved. Because of the aforementioned 
interplay between the scheme choice and the re- 
summation of running coupling singularity, this 
inconsistency is likely to affect strongly the small 
x behaviour. 

4.2. Coefficient function resummation 

The leading small x contributions to partonic 
cross sections are known to all orders for a small 
but increasing number of physical processes: they 
were first computed for heavy quark photo-and 
electroproduction in Ref. [21] (later extended to 
hadroproduction in Ref. [32]), they have been de- 
termined for deep-inelastic scattering in Ref. [27] 
and more recently for Higgs production 33] and 
the Drell-Yan processes [34] . 

Expressions for coefficient functions in the NLO 
of the double-leading expansion were already 
constructed in Ref. [7J, where, however, the run- 
ning coupling terms discussed in Sect. 12.31 were 
still left unresummed, thereby simplifying issues 
of scheme dependence, but at the cost of not re- 
producing the correct small x behaviour. How- 
ever, running coupling corrections to the resum- 
mation of coefficient functions also grow as £ — > 
co, analogously to the running coupling correc- 
tions to splitting functions Eq. they must be 
resummed lest physical observables develop un- 
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Figure 10. The matrix of coefficient functions 
with nt — 4, a s — 0.2. The curves correspond 
to NLO (dashed), NNLO (dot-dash), resummed 
(solid). The two different resummed curves (in 
the gluon sector) correspond to the MS (steeper 
at small x) and QoMS factorization schemes. 

physical singularities which leads to a spurious 
small x growth [50] . 

The resummation can be performed ex- 
actly [50] when the double-leading expression of 
the coefficient function is known in closed form. 
This is however not the case for the F2 deep- 
inelastic coefficient function in the MS scheme, 
for which only a series expansion generated by a 
recursion relation is available [27]. In such case, 
the dominant running coupling corrections to the 
coefficient function can be resummed through a 
divergent asymptotic expansion, which may be 
summed by the Borel method [17] . The ensu- 
ing resummed coefficient functions are displayed 
in Fig. [ITJ] Resummed coefficient functions were 
also presented in Ref. 20J (also including running 
coupling resummation) and are qualitatively sim- 
ilar, though a detailed comparison is hampered by 
the fact that the factorization scheme used there 
is different (DIS instead of MS). 

4.3. Parton distributions and structure 
functions 

Combining the ingredients discussed so far it is 
possible to determine resummed predictions for 
deep inelastic structure functions. Eventually, 
these should be used, together with resummed 




Figure 11. The ratio of the resummed (solid) or 
NNLO (dot-dashed) to NLO singlet quark and 
gluon distributions as a function of x at the scale 
Qo = 5 GeV (top) and as a function of Q at fixed 
x = 10~ 2 , 10~ 4 , x = 10~ 6 (bottom; the small- 
est x is the lowest curve in the resummed case 
and the highest at NNLO). The ratios are deter- 
mined assuming that the structure functions F2 
and Fl are kept fixed at the scale Qo = 5 GeV. 
The two different resummed curves correspond to 
the MS (smaller at small x) and QoMS factoriza- 
tion schemes. 

expressions for other physical processes, for the 
determination of parton distributions at the re- 
summed level. 

An estimate of the impact of resummation on 
parton distributions can be obtained by first com- 
puting the structure functions F% and Fl with 
some typical "toy" set of NLO parton distribu- 
tions (PDFs), and then assuming that the struc- 
ture functions are kept fixed at some scale: this is 
then enough to determine the resummed singlet 
quark and gluon distribution at that scale. The 
effect on PDFs is close to that which would be 
obtained if PDFs were determined from a fit to 
DIS data mostly clustered around that scale. Re- 
sults for the typical HERA (compare Fig. [J) scale 
choice of Q 2 = 25 GeV 2 are shown in Fig. [11] 
where we display the singlet quark and gluon dis- 
tributions as a function of x at this starting scale, 
and then as a function of Q 2 for various x values. 
Results are shown as a ratio of the resummed or 
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NNLO result to the NLO one. The impact of the 
resummation at the "HERA scale" is compara- 
ble to that of NNLO corrections, but it goes in 
the opposite direction: it tends to suppress the 
starting PDFs while the NNLO tends to enhance 
them. When evolving up the differences tend to 
be washed out because of asymptotic freedom. 

The PDFs displayed in Fig. [IT] can then be 
used to determine if-factors for resummed struc- 
ture functions: these were computed in Ref. |17j . 
There the factorization and renormalization scale 
dependence of the of results were also studied, 
and shown to be small. Indeed, once all ingre- 
dients are included, the resummed perturbative 
expansion of physical observables turns out to be 
as good as the standard fixed-order expansion. 

The dominant qualitative feature of the J^tT — 
factors is that resummation leads to a suppres- 
sion of the structure functions F2 and Fl at small 
x. The if-factors were compared to those deter- 
mined using the TW approach [20] in Ref. [21] . 
The main differences are related to the fact that 
TW if -factors can differ sizably from one even 
at x < 10~ 2 where the ABF result matches 
smoothly to the NLO one, and that the TW K- 
factors still show a marked scale dependence at 
large Q 2 ^ 10 4 GeV 2 where the scale dependence 
of the ABF result flattens out completely. These 
differences are likely to be due to the various fea- 
tures of the TW approach mentioned above: in 
particular, the lack of matching terms to large x, 
and the fact that the factorization scheme is not 
treated in a fully consistent way. Indeed, there 
appears [17] to be a significant cancellation of 
scheme dependence between evolved parton dis- 
tributions and coefficient functions, which is in- 
evitable spoiled if these are not determined in the 
same factorization scheme. 

if -factors can also be used to determine re- 
summed predictions in an approximate way, by 
first, computing the K factors as discussed above, 
but using a realistic underlying set of PDFs, then 
verifying that the values of the K factors are rea- 
sonably independent of the choice of underlying 
PDFs, and finally using the if-factors to correct 
NLO predictions obtained from any PDF set. 

In Fig. [T2] we show the structure functions F2 
and Fl at the scale Q 2 — obtained applying 
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Figure 12. The deep-inelastic structure functions 
F 2 and F L at Q 2 = M$y, computed at NLO, 
NNLO and resummed level. The parton distri- 
butions are from the NLO NNPDF1.0 set in the 
NLO case, while at the NNLO and resummed 
level they are computed using if -factors as dis- 
cussed in the text. The PDF uncertainty band 
is shown for the NLO prediction, while only the 
central prediction is displayed in the other cases. 

cither resummed or NNLO if -factors to the NLO 
prediction from the NNPDF1.0 parton set [35] . 
Interestingly, the suppression due to the resum- 
mation in the small x & 10 -4 region is larger than 
the one-sigma band due to PDF uncertainties. 

These results show that resummation is rel- 
evant for deep-inelastic scattering at a high- 
energy hadron-electron collider such as a collider 
based on the LHC (LHeC 36 ). It will be inter- 
esting to study the impact on precision LHC pro- 
cesses such as W and Z production of resumma- 
tion corrections, both from their direct effect on 
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partonic cross sections, and their indirect effect 
due to their impact on the extraction of parton 
distributions. 

In summary, we have seen that the impact of 
resummation is as large as that of NNLO correc- 
tions in the HERA region. The theoretical frame- 
work on which resummed results are based is now 
on a similar footing as that of standard fixed or- 
der perturbation theory, whose reliability is thus 
extended also to kinematical regions, relevant for 
HERA and the LHC, where NNLO results be- 
come unstable. 
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